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Abstract 

Elliott Lieb's ice-type models opened up the whole field of solvable models 
in statistical mechanics. Here we discuss the "commuting transfer matrix" 
T, Q equations for these models, writing them in a more explicit and trans- 
parent notation that we believe offers new insights. The approach manifests 
the relationship between the six-vertex and chiral Potts models, and between 
the eight-vertex and Kashiwara-Miwa models. 
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1 INTRODUCTION 

Elliott Lieb used the Bethe ansatz to solve the two-dimensional ice, F and 
KDP models in 1967 - all typical cases of the general "six-vertex" model 
whose solution was then given by Sutherland. - |5] This work, together 
with Onsager's famous solution 6 of the square-lattice Ising model in 1944, 
formed the basis on which the field of two-dimensional models rapidly grew. 
It led the author to develop the "commuting transfer matrix" of tackling such 
problems, notably the eight-vertex and hard-hexagon models. 
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This method begins with two steps. The first is to treat the problem 
in sufficient generahty that it contains an arbitrary parameter (a complex 
number) v such that the transfer matrix T is a non-trivial function T{v) of v, 
and 

Tiu)T{v) = T{v)T{u) (1) 

for all values of u and v. Thus T{u) and T[v) commute. The variable v is 
known as the "spectral parameter": for the six and eight-vertex models (and 
many other planar models) it is the difference of two "rapidity" variables 
associated with the horizontal and vertical directions of the lattice. 

The second step is to construct (if only implicitly) another matrix function 
Q{u) which also commutes with T{v), i.e. 

Q{u)T{v)=T{v)Q{u) , (2) 

and satisfies a matrix functional relation that is linear in T and Q, and ho- 
mogeneous in Q. For the six- vertex model this has the structure 

T{v)Qiv) = (l){v-X)Qiv + 2X') + (l){v + X)Q{v-2X') , (3) 

(p{v) being a known scalar function and A' a "crossing parameter". 
A further step is to show that Q{u) and Q{v) commute: 

Q{u)Q{v) = Q{v)Q{u) , (4) 

for all u, V. 

The final step is to go to a representation in which T{v),Q{v) are diagonal 
matrices, for all v. Then @ becomes a scalar functional relation for each 
eigenvalue, which can be solved. 

For the six and eight-vertex models, all these steps have been carried out 
by the author [71IHlini and written up in chapters 9 and 10 of ^U]. However, 
the construction of Q{v) given therein is implicit rather than explicit. 

Further, in the early papers |H1 El the author was concerned to focus 
on cases where there was no reason to believe the matrix Q{v) was in general 
singular. Obviously if it is zero the equations 0-0 contain no information 
on T{v). Only if it is non-singular (i.e. has non-zero determinant) for general 
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values of v can one expect to obtain all the eigenvalues of T(v) from (0). 
For this reason the author initially focussed on the "roots of unity" cases, 
where one can write Q{v) more explicitly. Later he was able to remove this 
restriction, and the working in ^U] is a quite general solution of the zero-field 
six and eight-vertex models. 

Here we re-present this working for the six-vertex model, giving a more 
explicit expression for the matrix Q{v).^ The relations © ~ @ are derived 
as the set of equations (|58|) , where the skew parameter r is exhibited and (for 
L even) the other skew or field parameter s can be given the value 0. 

The formulae ()27|). H29|) for the transfer matrix Qr{v) are interesting in 
that it is almost that of a trivial one-dimensional model: the only thing that 
stops this happening is a factor x'^^""'') in the Boltzmann weight function (|27() . 
A very similar property occurs in the three-dimensional Zamolodchikov model 
and its extension. jlH IT^ 

In 1990 Bazhanov and Stroganov IT3l showed how the recently-discovered 
solvable chiral Potts model could be obtained from the six-vertex model via 
the matrix Q{v). They used the general Q matrix discussed in section 6 
(also including the fields we mention at the end of our conclusion).^ In this 
approach the Boltzmann weight of the chiral Potts model first appears as an 
auxiliary function that enters the derivation of . Here we shall observe such 
a property. 

We emphasize that our working is for all values of the crossing parameter 

A or A'. This means that the spins in our spin formulations of the six- vertex 

and Q models are free to take all integer values from — oo to -|-oo. We use this 

convenient "spin language " to derive the relations, but then in section 5 we 

indicate how one can transform back to the arrow language of Lieb so as to 

ensure that the row-to-row transfer matrices are finite-dimensional. 

^Basically our results are more explicit versions of equations (9.8.15) - (9.8.37) of |1(J| . 

^It is the T2{tq) matrix of [21, except that it is rotated through 90°. Bazhanov and Stroganov 

rotated the lattice and started with the column-to-column transfer matrix of the six-vertex model. 

This leads to the usual row-to-row transfer matrix of the chiral Potts model. Here we shall not 

be extending the calculation that far, so revert to using the row-to-row transfer for the six-vertex 

model. 
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Only in the conclusion do we address the "roots of unity" cases, when 
A' = irmr/N, x^^ = 1, {m,N integers) and the spin on a given site takes 
just states. These cases are of course of interest. Fabricius and McCoy 
[T^ - ^Hl have studied these cases and have emphasized that the eigenvalues 
of T{v) are then degenerate.^ This provides a motivation for establishing 
the Q, Q commutation relations directly: if the eigenvalues of T^{v) are non- 
degenerate, the Q, Q relations are implied by those for T, T and T, Q. Further, 
it is this step that first introduces the weight function of the chiral Potts model. 
Only for the A^-state root-of-unity case does this model seem to be properly 
defined. 

Finally, we briefly indicate how the same approach applied to the eight- 
vertex model leads quite directly to the Kashiwara-Miwa model, As Hasegawa 
and Yamada have shown (Ml, the Kashiwara-Miwa model is a "descendant" of 
the zero-field eight-vertex model in the same way that the chiral Potts model 
is a "descendant" of the six-vertex model in external fields. 

This answers a question that has long puzzled the author: can one obtain 
a generalization of chiral Potts by using the Bazhanov-Stroganov method, but 
starting with the eight-vertex model instead of the six-vertex? The answer 
appears to be no. The model from which Bazhanov and Stroganov start is 
the six-vertex model in a field. Hasegawa and Yamada start from the zero- 
field eight-vertex model. There is no known solvable model that continuously 
generalizes and includes both these vertex models, so there is no reason to 
suppose that there exists a model which similarly generalizes and includes the 
chiral Potts and Kashiwara-Miwa models. 

It is true that the two vertex models intersect in the zero-field six-vertex 
model. The "descendant" of this is the Fateev-Zamodochikov model pij. 
which is indeed the intersection of the chiral Potts and Kashiwara-Miwa mod- 
els. 

The relationships between these six models are shown in Figure ^ 
^Fabricius and McCoy use the word "incomplete" in a way quite different from its usual sense 
in mathematical physics ^1 1213 1^ ■ 
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2 TRANSFER MATRICES AND THE STAR- 



TRIANGLE RELATION 

Here we make some general observations to which we shah refer in subsequent 
sections. 

We shah consider various "interactions-round-a-face" (IRF) models on the 
square lattice £ of TV sites. In each model, each site of the lattice carries a 
"spin" . Let a, b, c, d be the four spins round a face, arranged as in Figure 12 
and let W{a, b, c, d) be the Boltzmann weight of this spin configuration. Then 
the partition function is 

^ = En^(«'^'C,d) , (5) 

the product being over all faces and the sum over all allowed values of all 
the TV spins. ^ There are usually restrictions on the spin values: here we shall 
require that on some (not necessarily all) edges the two adjacent spins 
(Tj, aj differ by unity: 

aj = ai±l . (6) 

In particular, we require that ^ be satisfied on all horizontal edges. Let 
the lattice have L columns, i.e. L faces per row, and let the spins in any 
particular row be ai,a2, ■ ■ ■ , (Tl+i- Then we impose skewed cyclic boundary 
conditions: 

(TL+i = cTi + r , (7) 

where r is the "skew parameter" - an integer with value between —L and L 
such that L — r is even. It may vary from row to row. 
We shall also use the spin differences: 

Oj = o-j+i - (^i ■ (8) 

we see that each Ui takes the values +1 or —1, and 

ai H h az, = r . (9) 

'^Some of the spins may need to be fixed to avoid this sum being infinite. For the six-vertex 
model just one spin on the lattice needs to be fixed - say the one at the bottom-right corner. For 
the "Q" model we shall discuss we need to fix one spin in every row - say the left-most spin. 
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Let a = {o"!, . . . , cTi+i} be the set of spins in one row, and let a' = 
{a[, . . . , cr'^^i} be the spins in the row above, as in Figure |21 Then the row- 
to-row transfer matrix T is the matrix with elements 

T^y = W{ai,a2,(T2,a[)W{a2,a3,a3,a2) . . .W{aL,crL+i,a'L^i,a'L) ■ (10) 

If the lattice has M rows, with cyclic top-to-bottom boundary conditions, it 
follows that 

Z = Trace . (11) 

Obviously T depends on the skew parameter r. If r is the same for both 
rows, we may write T as . If it is different, so that Ci+i = ai + r and 
(T^_j_]^ = a'l + s, we may write T as T^*. 

Commutation 

Now consider two models, with different weight functions Wi and W2- Let 
their transfer matrices be Ti, T2 (for the moment we suppress the superfixes 
r,s). The product T1T2 has entries 

[TiT2]a,b = Trace MiM2---Ml . (12) 

Here Mj is a a matrix which depends on four outer spins Oj, aj+i, bj, 6j+i and 
has entries 

[Mj]c,d = Wi{aj,aj+i,d,c)W2{c,d,bj+i,bj) . (13) 

The entries of T2T1 are similar, but Mj is replaced by Mj, in which Wi and 
W2 are interchanged. For Ti to commute with T2, i.e. for T1T2 = T2T1, we 
need there to exist invertible matrices Pi, ... , Pl+i such that Pj depends on 
the outer spins aj,bj (and no others), and 

M; = Pf'MjP,+i , j = l,...,L . (14) 

Remembering that the matrix Pj depends on aj,bj, write its entry (c, d) 
as W3{c,aj,d,bj). Then, rewriting (fTl|) as MjPj^i = PjMj, we obtain this 
condition explicitly as 

Wi {b, c, g,a)W2 (a, g, e, /) (5, c, d, e) = 

9 
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Wsia, b, g, f)W2{b, c, d, g)Wi {g, d, e, /) (15) 

9 

for all allowed values of the external spins a, b, c, d, e, f. 

For interaction-round-a-face models, (|15|) is the "star-triangle" relation. It 
is depicted graphically in Figure 21 

For cyclic boundary conditions, Pl+i = Pi and H15|) . together with the in- 
vertibility of Pj, is sufficient to ensure that T1T2 = T2T1, i.e. Ti, T2 commute. 

However, for our more general skewed boundary conditions we do have a 
problem. Including the skew parameters, the commutation relation in general 
becomes 

T^^T^' = TfTl'' . (16) 

Here r is the skew parameter of the lowest row, s of the uppermost, and t, t' 
of the intervening rows (the ones whose spins are summed over in the matrix 
multiplication). Note that t is not necessarily the same as t' . If Pi has entries 
Ws{c, ai, d, 61), then Pl+i has entries Wsic + t,ai + r,d + t' ,bi + s). If these 
entries are equal then (|16|) holds true. 

In fact they are not always equal: in section 4 we encounter a case where 

W3{c + t,a + r,d + t',b + s) = x'^""-"^ W3ic,a,d,b) , (17) 

X being the fixed parameter defined in pij). The effect of this is to insert an 
extra factor x^^'^'^~'^'^^ into the sum over ci,. . . ,cl that is the matrix product 
T1T2. This is equivalent to replacing Ti therein by D^^TiD'^, where D is a a 
diagonal matrix with entries 

Day = x'''5{ai,a[)---6{aL+ua'L+i) . (18) 

(Note that here a and a' necessarily have the same value of the skew parameter 

<7L+1 - o"i.) 

Thus instead of ((TB|) we obtain 

D-'Tl^D'Tl' = TI^'t(' . (19) 
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3 THE SIX- VERTEX MODEL 



As considered by Lieb and Sutherland, the six-vertex model is a model on 
the square lattice where one puts an arrow on every edge. At each site one is 
required to satisfy the "ice rule" that there be two arrows pointing in and two 
pointing out of the site. There are then six possible configurations of arrows 
at a site, as indicated in Figure |SJ 

For our present purposes, we wish to express this in interaction-round-a- 
face language. We do this by putting spins on the faces of the lattice so that 
as one goes from one face to a neighbouring face, if the intervening arrow 
points to one's left (right), then the spin value increases (decreases) by unity. 

As one walks around a site, the ice rule ensures that there are two arrows 
pointing to the left, and two to the right, so one does indeed return to the 
original spin value. 

It follows that a spin configuration is allowed iff every pair of adjacent 
spins (horizontal and vertical) differs by unity, i.e. satisfies ©. If one of the 
spins is fixed (say the one at the lower-right-hand corner) , then there is a one- 
to-one correspondence between allowed arrow and allowed spin configurations 
on the lattice. 

Finally, we go to the dual of the arrow lattice. The spins now live on sites 
(rather than faces), and Figure El becomes Figure El 

The spin differences mentioned above now describe the vertical arrows 
in a row: Oj = +1 if the arrow between spins ai and Cj+i is pointing upwards, 

= —1 if it is pointing downwards. We identify the horizontal arrow between 
(Ti and a'l with the arrow between cjl+i and (^•^- ™ arrow formulation 
we use the usual cyclic boundary conditions). Then the ice rule implies that 
if © is satisfied for one row, then it is satisfied for all rows, with the same 
value of r. This is the "conservation of up and down arrows" property of the 
six- vertex model. ^ - |Sj 

Let wi, . . . ,WQhe the Boltzmann weights of the six configurations in Figure 
Eland FigureEl Here we consider only the "zero-field" six-vertex model, which 
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is invariant under reversal of all arrows, so we take: 

Wi = W2 , W3 =W4 , W5 = Wq . (20) 

The overall normalization of wi,W2, 103 plays a trivial role in the calcu- 
lations, so without loss of generality we can define two parameters A, v so 
that 

W2 = sinh ^{X — v) , W4 = sinh ^{X + v) , wq = sinh A (21) 

as in (9.2.3) of 10 . It follows that the Boltzmann weight function of the 
zero-field six vertex model is 

W{a,b,c,d) = WQ^{v\a,b,c,d) = 5(a, c) e(^+'')(''+'^-''-'=)/^ sinhi(A-?;)+ 

5(6,d)e(^-^)(''+"-''-'^)/^ sinhi(A + i;) . (22) 

We regard A (the "crossing parameter") as a given constant and v as a 
variable (the "spectral parameter"). We therefore write the six- vertex model 
transfer matrix as T{v) = T^{v) and define it as in (|lfl|) : 

[^''(^)]<T,<7' = WQy{v\ai,a2,a'2,a[)...WQy{v\aL,(TL+i,a'L+i,(7L) ■ (23) 

The lower and upper spin sets a, a' have the same value of the skew parameter 
r. 

The star-triangle relation ()15j) is satisfied if 

Wi = W6,{v) , ^^2 = 1^6.(^^0 , W3 = We^iv' - V - X) 

for all values of v and v'. Also, taking r, s, s' , t = r, we see that (|17|) is satisfied 
by (EU, with x = l. 

The six-vertex model transfer matrices therefore commute: 

T'-{v)T'-{v') = T''{v')T''{v) . (24) 

Multiplying WQv{a,h,c^d) by a field factor ^'^^^ is equivalent to post- 
multiplying T'^[v) by a diagonal matrix with entries 

Sl,a' = Cicy) n '^(^J'^i) ' (25) 
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where 

This in turn corresponds to introducing an electric field acting on the upper 
vertical arrows, or equivalently on the upper pairs of horizontally adjacent 
spins. Clearly commutes with T^'{v): 

S'T'iv) = T''{v)S' . (26) 

4 THE "Q" MODEL 

The step taken by Bazhanov and Stroganov JH] was to look for another 
solution of ()15|). in which W3 remains unchanged, but Wi, W2 are altered so 
that they are no longer weight functions for the six-vertex model, but for a 
new "Q" -model. These functions can in fact be obtained (after changing the 
normalization) from (9.8.15) - (9.8.23) of ^U]. They are 

WQ{v\a,b,c,d) = ex.p[X'd{a-b)/2-{X + v){b-a){c-d)/A] , (27) 

where 

A' = A-ivr , (28) 

as in eqn. 9.3.7 of ^Uj. 

Horizontally adjacent spins must still satisfy ©, so |6 — a| = |c — d| = 1, 
but now there are no restrictions on the differences d — a or c — b of vertically 
adjacent spins. Hence the skew parameters r, s of the lower and upper rows 
in Figure 01 can be different. 

We define the corresponding transfer matrix Q^^{v) as in (fT(l|l : 

[QR{v)Ua' = WQ{v\ai,a2,(T2,a[) ■ ■ ■WQ{v\aL,crL+i,(^'L+i,(^'L) ■ (29) 

T, Q commutation. 

In the star-triangle relation (|15|) we now substitute 

Wi = Weviv) , W2 = Wq{v') , W3 = Wq{v' - V - X) . (30) 
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To be consistent with the definitions of Wq^ and Wq, we require that spins 
Unked directly by edges in Figure ^ differ by one, except for the edges {a,f), 
{g,e), {c,d) on the Ihs, and edges (a,/), {b,g), {c,d) on the rhs. This means 
that the six external spins break up into two sets: (a, b, c) and (d, e, /). One 
is free to independently increment all the spins in either set by unity. We find 
that (Uni) is satisfied by 

For the T and Q/j, matrices to commute, we also need to check the the 
auxiliary condition (|17j) . Since T necessarily relates two rows of spins with 
the same skew parameter, we must have t = r,t' = s therein. We find that 
H17|l is satisfied, but x therein is not unity: instead 

x = e^'/2 = _ieA/2 _ (31) 

We therefore obtain the modified commutation relation with Ti re- 
placed by T^{v) and T2 by Q^(u'), i.e. 

D-^T\v)D^Q'^{v') = Qn{v')T\v) . (32) 

Here v and v' are arbitrary, D is defined by (|18() and ()31() . The notation is 
threatening to become confusing: note that r and s here are merely superfixes, 
except that is actually the diagonal matrix D raised to the power s. 

The r, Q functional relation. 

Our working in section 2 depended on the matrices Pj being invertible. This 
will cease to be so if VF3(a, 6, c, d) factors into a function independent of a 
times a function independent of c. For the derivation of (|32|) this happens 
when v' = V. 

This does not mean that the commutation relation fails - it follows by 
taking the limit v' ^ v. It does mean that there is then additional information 
in the relation H15() . The Ihs then involves d only via the simple factor x^^^~^\ 
while the rhs involves a only via x^^^-^^K 

Going back to equations (fT2|) - (|TI|) . define vectors ^^i, . . . with elements 

= X^^^'^-"^) . (33) 
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Then from the above observations about the star-triangle relation 

Mj^j+i = hjCj , (34) 

where 

hj = h{aj,aj^i,bj+i,bj) 

is a scalar factor dependent (as are Mj and on the lower and upper spins 
Qj , Oj+i , bj , . In fact 

h{a,b,c,d) = J2WUv\a,b,g,f)WQiv\f,g,c,d)x<3~'^^-''^f~-^ . (35) 

9 

Here W^y, Wq are defined by ((2H), (HTJ. The rhs is necessarily independent 
of the extra spin /. 

A direct calculation of (|35|) reveals that 

h{a,b,c,d) = smli^{\-v)WQ{v + 2X'\a,b,c,d) (36) 

Since Wi in ()13() is the six-vertex model weight function, c is aj it 1 and 
d is Oj+i ± 1. Hence each Mj is a two- by-two matrix. Now take Pj to be 
a two- by-two matrix whose first column is ^j. The second column can be 
chosen arbitrarily, so long as the choice makes Pj invertible. It is convenient 
to choose det(P;) = 1. 

Then H34|) ensures that 

Pf'M,Pj+i = Mj , (37) 

where Mj is an upper-triangular two-by-two matrix with upper-left entry 
hj. Equating determinants, we find that the lower-right entry of Mj is 
h' {aj, aj^i,bj^i,bj), where 

h'{a,b,c,d) = smli^{X + v)WQ{v-2X'\a,b,c,d) . (38) 

Again, we have to worry about our skewed boundary conditions. In gen- 
eral we do not have Pl+i = -Pij but rather Pl+i = T^Pi, where P is a diagonal 
two-by-two matrix with entries 3;*(^~"i^(5(c, d). As in we must therefore 
replace Ti by D^^TiD'^, which is equivalent to replacing the rhs of (|12|) by 
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Trace VMi ■ ■ ■ . Then the similarity transformation (|37|1 reduces the ma- 
trix product to upper-triangular form, so the trace becomes the sum of two 
products, the first over all upper-left elements of the Mj, the second over all 
lower-right. We obtain 

D-'T\v)D'Q''^{v) = <P{X-v)Q"]^{v+2X') + (j){X+v)Q]^{v-2X') , (39) 

where 

cl){v) = [sm{v)]^ . (40) 

Q, Q commutation. 

Now we consider the conditions for two transfer matrices Q^^{v), with different 
values of v and possibly different values of r, s, to commute. This is the step 
at which one first encounters the chiral Potts model. 

As we have seen, the six-vertex model transfer matrices T satisfy the com- 
mutation relation (|24|). and the Q matrices satisfy the commutation relation 
(|32j) with T. At first sight this would appear to imply a corresponding commu- 
tation relation between the Q's themselves. If the T have unique eigenvalues, 
this is so. 

However, Fabricius and McCoy have emphasized that for the special case 
when our parameter rr is a complex root of unity, then the eigenvalues of the 
T matrices are degenerate, and interesting algebraic structures emerge.^ It is 
therefore desirable to establish the Q, Q commutation relations directly. 

We immediately strike a difficulty with our spin language. Up to now we 
have only considered matrix products that involve at least one six-vertex T 
matrix, e.g TQji. Each element of this product corresponds to two rows of 
faces of the lattice, the lower with transfer matrix T, the upper with Q^. 
There are three rows of spins on sites, those in the lowest and uppermost are 
given, the ones in the middle row are to be summed over. Since at least one 
of the matrices is T, the ice rule (all horizontally and vertically adjacent spins 
^The eigenvectors of T arc then of course not unique, so it is not surprising that (|39|) then 
defines the eigenvalues of T uniquely, but not those of Q - a fact that Fabricius and McCoy seem 
to find remarkable. 
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differ by one) ensures that that there are a finite number of terms in this 
summation. 

However, if both matrices are Qr matrices, one can choose one of the 
spins in the middle row arbitrarily. At this stage we are allowing spins to take 
all integer values, so there are an infinite number of choices and an infinite 
number of terms in the summation. The matrix product is not defined. 

To overcome this, we anticipate the next section, where we go from spin 
language back to "horizontal spin difference" or "vertical arrow" language. 
If we can arrange a QQ product so that the terms in the sum it represents 
are unchanged by incrementing all the spins in the middle row by unity, then 
we can write this sum as a sum over only the differences of such spins. This 
is equivalent to arbitrarily fixing one of the spins in the middle row, say the 
left-hand spin. 

From (|27)) and (|29j) . incrementing by unity all the spins in the lower row 
of the matrix Q^j leaves the elements of that matrix unchanged. It is not so 
the upper row: that introduces a factor x"*" coming from the exp[A'(i(a — 6)/2] 
term in (|27|) . We can overcome this by instead using the matrix 

Q'^liv) = x'^D'Q'-^{v)D^ (41) 

which has an extra factor in H29|) . This is in turn equivalent to 

multiplying ((2Z|) by 2;'"=-"'^, i.e. to replacing X'd{a-b)/2 by X'b{c-d)/2. Thus 

[Ql iv)Ua' = WQ{v\ai,a2,(J2,a[) ■ ■ ■WQ{v\aL,(TL+i,(T'L+i,crL) , (42) 

where 

WQ{v\a,b,c,d) = exp[X'b{c-d)/2-{X + v)(b-a){c-d)/A\ . (43) 

This matrix Ql corresponds to the Ql of (9.8.27) of Ref. but with the 
normalization changed and both A and v negated. 

This matrix Q^j^ is unchanged by incrementing all the upper spins by unity, 
and we can define the product 

r'%v,v') = Q'^iv)Q'^{v') (44) 

as above. 
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Let {bi, . . . , be the spins that are summed over in this matrix prod- 

uct, i.e. the spins on the middle row, above Q^j^iv) and below Q^^{v'). They 
enter only via their differences f5j = fej+i — bj. The sum is over all values (±1) 
of the subject to the skew condition 

A + ■ • • + /3l = t . 
We can remove this restriction by forming the generating function matrix 

F'-\z\v,v') = Y,z'r''{v,v') (45) 

t 

the sum being over all values of t from —L to L, with L — t even. Then 
F'^^{z\v,v') is given by the unrestricted sum over the Pj. 

Being a transfer matrix product, the elements of F^^{z\v,v') are given in 
the first instance by an expression of type dJ). However, now the elements 
[-^j]c,(i of Mj depend on the two intermediate spins c, d only via their difference 
m = d — c = ±1. The trace in (|12|) simplifies to a sum over all the L spin 
differences. These are now independent, so we can perform the m summation 
individually for each Mj. The rhs of H12|) then becomes a product, as in 
Hence the elements of F^'^{z\v,v') are given by the rhs of with W 
replaced by Wp, and 

WF{v,v'\a,b,c,d) =Y,z"'WQ{v\a,b,f + mJ)WQ{v'\fJ + m,c,d) , (46) 

m 

the sum being over m = —1 and +1. The rhs is independent of / and from 
([27|l, (pl) 

Wf{v, v'\a, b, C,d)=Y^ ^m^mib-d) ^m[{X+v)m{a-b)+{X+v'){d-c)]/A _ (47) 
m 

Here v, v' are independent variables. 

We look for an auxiliary function P{v,v'\a,b) such that 

WFiv',v\a,b,c,d) = P{v,v'\d - a)WFiv,v'\a,b,c,d)/P{v,v'\c- b) , (48) 

since then (apart possibly from boundary conditions) the rhs of (|lfl|) is unal- 
tered by interchanging v with v', as the P factors cancel out of the product 
on the rhs. This in turn implies a commutation relation between Q^i^{v) and 

QtW). 
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There are four cases to consider in (|48|) . b — a = ±1 and c — d = ±1. If 
b — a = c — d, then ()48() is automatically satisfied. The other two cases yield 
just one condition on the function P, which can be written 

Piv,v'\j + 1) ^ ze^^ + x^^e^^' 

P{v,v'\j-l) ze2^' +x2ie2'^ ' ^ ^ 

Hence P{j/2) has the same structure as the chiral Potts weight function 
Wpg{j) in eqn (2) of [^S], with to = x'^. However, we emphasize that at 
this stage we are not requiring to be a root of unity, so the rhs of H49|) is 
not necessarily periodic in j. 

Substituting the form for W = Wp into (fTU]) . the P- functions can- 
cel except for a boundary factor P{v,v'\a'^ — ai) / P{v^v'\a'j^_^^ — aL+i)- Re- 
membering that cJi+i = (Ti + r and o"^_|_x = + this has the form 
P{v,v'\d — a) / P{v^v'\d — a + s — r). This is unity if s = r, so we have es- 
tablished the symmetry property F^'^ {z\v' ^v) = F'^^{z\v,v'). This is true for 
all z, so the coefficients P^^'{v,v') in the Laurent expansion of F"'(z\v',v) are 
also symmetric in v,v', for all t. Hence from (|44|) 

Ql{v)Q%{v') = Ql{v')Q%{v) ■ 
Prom (|41|) . replacing i by s, this becomes 

Q'j^{v)D^Q'^{v') = Q'i{v')D^Ql{v) . (50) 



5 SUMMARY OF THE RELATIONS 

Neither Q^^{v) nor Q'j^{v) has the desirable property that it is unchanged by 
incrementing all the spins in either the lower or the upper row in Pigure |31 
by unity. Prom the discussion preceding ()41() . the matrix that does have this 
property is 

Q-%v) = Qr{v)D'- = x-^^'D-'Q^iv) . (51) 

It is convenient to introduce a generalization of the six- vertex model trans- 
fer matrix T'"(u): 

T^{v) = D-'T'-{v)D' . (52) 
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This is the transfer matrix of the six-vertex model with skew parameter r 
and a field of weight x'^^~""^ acting on the two vertically adjacent spins (or 
equivalently the horizontal arrow between them) at the left-hand end of the 
row. Then the commutation and functional relations (|24() . H32|) . ()39() . (jSUj) for 
the T and Q matrices become 

T:{v)T:{v') = Tl{v')Tl{v) , (53) 
T:{v)Q^\v') = Q'^v'mv) , (54) 

T:iv)Q'^%v) = <PiX-v)Q'-'{v+2X') + cj){X+v)Q''{v-2\') , (55) 

Q''{v)Q''iv') = Q''{v')Q"'{v) (56) 

for all values of the skew parameters r, s. Note that nowhere in this paper do 
we use any implicit summation convention over skew parameters. 

Square matrix form of the relations 

We have uncovered quite an interesting structure in the T, Q relations for the 
six-vertex model, and we could continue to work with the above relations (|53|) 
- (|56|) . However, they are still not quite in the form of the T,Q (or V,Q) 
relations given in (9.2.5), (9.8.40), (9.8.42) of Ref. [TOl. These involve the 
zero-field six- vertex matrix T^{v) and a square invertible matrix Q{v), for 
given r. 

To derive these square matrix relations, define, using (jSB)) . 

Qliv) = Q'\v) Q-(t;o) = Q'^v,) Q^'iv) (57) 

vq being some fixed value of v. 

Then from (|53j) - (|56|) we readily obtain the desired final square matrix 
relations 

t:{v)t:{v') = n{v')T:{v) , 
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TI{v)Q:{v') = Ql{v')Tl{v) , 



(58) 



T:iv)Qliv) = cPiX-v)Q:{v+2X') + cl>{X+v)Q:iv-2X') 



Ql{v)Ql{v') = Ql{v')Ql{v) 



true for all the allowed values — L, 2 — L,4 — L, . . . , L of r and s. 

From our remarks in the next sub-section, we expect Ql{v) to be invertible 
provided Cs> Cr- For an even number L of columns, this will always be so if 
we choose s = 0. Then Tq{v) = T{v) is the usual zero-field six- vertex model 
transfer matrix and Qq{v) = Q{v) is the associated Q-matrix used in - ^01 
. We have derived the relations O - ®) except that they can be and usually 
are thought of in the "arrow language" of the next sub-section. 

Transformation to arrow language 

We have used the language of spins on sites to derive these relations, but the 
r, Q matrices can be thought of conveniently in the alternative spin-difference 
or "arrow" language, transforming from the aj to the aj by ©. The equations 
(ITHl) . (ITtI) . (PH) . ((SU) define the elements {a, a') of Q'''. They depend only on 
{a, a'), so we take the elements to be labelled by a, a'. Since the aj satisfy 
©, and the a'j satisfy © with r replaced by s, the matrix Q^^ is of dimension 
Cr by Cs, where 



and Ql{v) is a square matrix, of dimension Cr- 

From this point of view the six-vertex model transfer matrix is slightly 
more subtle, since we are not free to increment all the spins in the upper row of 
Figure El independently of those in the lower row. Only the face configurations 
of Figure are allowed. If the spins in the lower row are known , then there 
two choices for those in the upper, depending on whether a'l = cJi + 1 or 
a'l = (Ti — 1. However, the above matrix products involving T always involve 
a sum over these two choices. The end result is that the above relations ()53|) - 
are also true in the arrow language provided we sum over the two choices 
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of a[. From (UHl), (gSI), dSH) we obtain 

(59) 

the sum being over (5 = ±1, where (3 = a'l — ai. The rhs is uniquely defined 
by ai,...a^. This is the usual six- vertex model transfer matrix in arrow 
language, with an extra weight x'^ {x"") if the left-hand horizontal arrow 
points to the left (right). It is of dimension by C^. 

Transformation to all Tl{v) diagonal 

We assume that the matrices TJ(?;), Ql{v) are diagonalizable. Then from 
the set of relations H58|) . there exists a non-singular matrix VI, independent 
of V, such that the transformation TJ(w) — s- V^T^ {v){Vl}~^ reduces both 
TJ(f) and Ql{v) to diagonal form. If we also transform Q'^'^{v) by Q^'^{v) 
VIQ^'^ {v){V^}~^ , then the relations - and ((^ are unchanged. 

We do not have a proof for general L, but for L = 2 to L = 6 we have 
verified with Mathematica that the matrix product Q^^{v)Q^^{v') is invertible 
if Cs > Cr, for all such values of r and s and for arbitrarily chosen values of 
V, v' , X. This implies that the rectangular matrix Q^^{v) is then of "full" rank, 
having the maximum rank possible for a matrix of its size: 

rank of Q""' = min(a,C,) . (60) 

From the transformed form of ()54() this implies that if Cs > C^, then 
the eigenvalues of TJ(t>) are all contained in the eigenvalues of T^{v), which 
is an interesting result of which the author was previously not fully aware. 
If the eigenvalues of T^{v) are non-degenerate, (|54() further implies that one 
can order the eigenvalues so that Q'^^{v) transforms to a diagonal (but not 
necessarily square) matrix. Even if they are degenerate, then the fact that 
the transform of Ql{v) is also diagonal may imply this property of Q'^^{v). 
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Structure of the eigenvalues 

li Cs >= Cr, then it follows from ()27p and H29|) that for r >= all elements 
of Q'^{v) are of the form 

^v{r-s-L)/4 X { polynomial in e" of degree {L-r)/2} ; 
while for r <= they are of the form 

gi;(s-r-L)/4 ^ | polynomial in e^ of degree (L + r)/2 } . 

li Cs <= Cr, then the above statements are true if one interchanges r with s 
therein. (All four statements are actually true for all r and s, provided one 
allows some of the initial and/or final coefficients of the polynomial to be zero, 
which means that some of the vj below will be infinite.) 

From and (|F7|) . the same must be true for the elements of Ql{v). 
Also, since the eigenvector matrices "PJ are independent of v, it is true of 
the diagonalized form of Ql{v), i.e. of its eigenvalues. If Cg >= Cr we can 
therefore write any eigenvalue as 

(L-r)/2 

if r >= 0, or 

(L+r)/2 

Ce"'/^ n sinhi(t;-?;j) 

if r <= 0. Here C, vi,V2, ■ ■ ■ are constants, independent of v. 

Similarly, the eigenvalues of Tg(v) necessarily have the structure 

L 

C' Y[ sinh ^{v — Wj) , 

i=i 

C, wi,W2, ■ ■ ■ being other constants. 

We can substitute these forms into ()58|) and in principle solve for the 
various constants. Setting v = Vj gives the usual Bethe ansatz equations. 
Alternatively, one can expand the functions as polynomials in e^ and equate 
coefficients. 
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6 A MORE GENERAL FORM FOR Q 

We remark that (|47|) is the Boltzmann weight of the "r2" model rotated 
through 90° - eqns. (3.44), (3.48) of Ref. iMj. If instead of ^ we take 

Wi = We^iv) , Wa = Wf{v', v") , W3 = Wf{v' -v - X,v" -v- X) , 

then the star triangle relation ()15|) remains satisfied, for all v,v',v". 

This is the quite general solution of (|15|) found in 1990 by Bazhanov and 
Stroganoff.|13j There are further trivial factors that can be introduced into the 
Boltzmann weight Wpiv, v'\a, b, c, d) without affecting {TSJ, notably fi^^"" fi'^ '^ , 
where fii,^2 are arbitrary constants (the same for both W2 and W3). This 
factor merely multiplies the transfer matrix F^^{z\v,v') by ^1/^2- 

From (grj), (j331), (jlSI), (EU and the discussion before (fSU)) . 

F''{v,v') = F'\v',v) = F'''-{z\v,v') =Y,z^x'*D*Q''\v)Q*''{v')D~* . 

t 

It does indeed follow from that 

T'-{v)F''{v',v") = F'-{v',v")T'^{v) , 

for all v,v' ,v" . This is in agreement with the fact that Wpiv' , v") satisfies ()15(l . 
We stress that this is a "spin language" result: F'^[v\v") is not unchanged 
by incrementing all the spins in either the lower or upper row of Figure El by 
unity. 

7 CONCLUSION 

We have obtained the "T, Q" functional matrix relations for the zero-field six- 
vertex model (and for other special values of the horizontal field) , presenting 
the working in a way that we hope is more explicit and transparent than those 
used previously. As a check on our reasoning, we have verified the relations 
H53|) - ()56() explicitly on a computer, using Mathematica, for L = 2 to L = 6. 

In some ways this method appears to be less general than Lieb and Suther- 
land's original Bethe ansatz solutions: it is not obvious how to extend it to 
the more general situation when there is an arbitrary horizontal electric field 
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(which corresponds to a vertical spin field), or even to solve the zero- field 
problem with an odd number of columns. Still, it should be remembered that 
the commuting transfer matrix method was originally developed as a means 
of solving the eight-vertex model, for which there was then no Bethe ansatz. 

We emphasize that all the above working is for the general case, when 
there are no restrictions on the real or complex "crossing parameter" A. 

The "root of unity" cases. 

A very interesting special case occurs when x*^ is an A^th root of unity, i.e. 

X' = im-jr/N , x^^ = 1 , (61) 

m, N being integers and k being odd. These cases have been studied by 
Fabricius and McCoy, and they have shown that the eigenvalues of T^{v) are 
then degenerate. - jTH] 

One can then regard the lattice spins cij, a'^ as restricted to N values. For a 
given site these will either be odd or even, e.g. 1,3,..., 2N — 1. Horizontally 
adjacent spins differ by ±1 modulo 2N , so we continue to take the aj in 
((SI) to be strictly ±1, but interpret the difference on the rhs to modulo 2N . 
Similarly, the skew parameter r, given by ((J)), is to be interpreted modulo 
2N . One can then repeat all the above working. (There may have to be some 
minor modifications, such as multiplying <5^(f ) by e"^^/^ so as to ensure that 
incrementing the spin d in (|27|) by 2N leaves the element of Q^^{v) unchanged. 
The horizontal spin differences a — b, c — d therein are to be kept as strictly 
±1.) 

There will be one significant difference: in the above we have specialized 
to the case when no horizontal spin fields (corresponding to vertical arrow or 
"electric" fields) are applied to the lower or upper rows of the T or Q matrices, 
for instance including a factor fi^-" in (|?7|) . We could have introduced such 
a factor as it does not affect the star-triangle relation (|15j) . (|3fl|) . It merely 
post-multiplies <5^(f) by the diagonal matrix 5^ defined in section 3. Since 
5^ commutes with TJ(u) and Qs{v), the equations would have been 
unchanged: no increase in generality would have been achieved. 
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For the root of unity cases this ceases to be true: fi^ is not in general 
the same as /i^'"*"^^. now commutes with TJ(f), but not with Ql{v). For 
generahty it is therefore important to include such factors in H27() . This leads 
to an extra constant term C on the rhs of H49|) . From ()48l ) . remembering that 
\a — b\ = \c — d\ = 1, the function P{v, v'\j) is defined only for either j always 
odd, or j always even. Taking j to be odd, from H49() . 

P{v,v'\2n-l) = Y[C \A , (62) 

where 

CJ = , OJ^ = 1 . 

The function P{v, v'\j) must now be periodic in j of period 2N , which will 
be so iff 

With a correct identification of the various parameters, we see that 
P{v,v'\2n — 1) is precisely the A^-state chiral Potts model weight function 
Wpq{n). Further, the factor C arising from the applied fields is important: 
without it one would have z = e^'^l^ and the model would reduce to the 
critical Fateev-Zamolodchikov model. [21] 

Kashiwara-Miwa model as a descendant of the eight- 
vertex model 

Just as the chiral Potts model can be regarded as a "descendant" of the six- 
vertex model in a vertical arrow field, so Hasegawa and Yamada have shown 
that the Kashiwara-Miwa model can be regarded as a "descendant" of the 
zero-field eight-vertex model. (23| I22j Note that neither vertex model includes 
the other: they intersect in the zero-field six-vertex model. Correspondingly, 
the chiral Potts and Kashiwara model intersect in the Fateev-Zamolodchikov 
model. 

In the same way that our present approach leads very directly from the 
six-vertex model to the chiral Potts model, so does it lead from the eight- 
vertex model to the Kashiwara-Miwa model. For the eight-vertex model, the 
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function Wp{u,v\a,b,c,d) is the multiplicand of eqn. (10.5.27) of 10 , the 
Sj,s'j,aj,aj therein being 

Sj = s + Xd , s'j = s' + Xa , aj = c — d , a'^ = h — a , (63) 

and s,s' are arbitrary parameters.^ 

Here we established the Q, Q commutation by showing that F'^'^(z\v, v') = 
F^'^ (z\v' , v) . Similarly, in ^] the author established the relation (10.5.29) by 
showing that the product (10.5.27) is a symmetric function of v and v'. That 
argument was by inductive reasoning on the entire product. If instead we look 
for a solution for the function P in H48() . we find that the P- factors depend 
on a, d (or 6, c) not just via their difference, but now P{v, v'\a, d) is a product 
of a function of d — a difference and a function ol a + d: 

P{v,v'\a,d) = Pi{v,v'\d - a)P2{v,v'\a + d) . 

Focussing on the root of unity case, when A' = A — 2ii^ = 2im K/N, m and 
being integers, and requiring that P{v, v'\a + 2N, d) = P{v, v'\a, d+ 2N) = 
P{v, v'\a, d), we find that P{v, v'\a, d) is the Boltzmann weight function of the 
Kashiwara-Miwa model. Since the eight-vertex model transfer matrix does 
not commute with S"^, it does not break up into diagonal blocks with different 
r and there is no analogue of the arbitrary factor C in (|62j) . 
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Figure 1: The three vertex models and their respective "descendants". The central 
models are special cases of the ones to their left and to their right. 
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Figure 2: Four spins round a face of the square lattice. 
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Figure 3: A typical row of L faces of the square lattice, with spins cJi, . . . , (Tl+i on 
the lower row of sites, cr^, . . . , cr^+i on the upper. 
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Figure 4: The generalized star-triangle relation. 
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1 2 3 4 5 6 

Figure 5: The six arrow configurations of tlie six- vertex model. 



a— 1 a a+1 a a—1 a— 2 a+1 a+2 a—1 a a+1 a 



a a+1 a a—1 a a—1 a a+1 a a—1 a a+1 



Figure 6: Tlie six configurations of tlie six- vertex model, expressed in terms of spins 
on the dual lattice. 
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